Entangled Spin-Orbital Phases in the d 9 Model 
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We investigate the phase diagrams of the spin-orbital d 9 Kugel-Khomskii model for a bilayer and 
a monolayer square lattice using Bethe-Peierls- Weiss method. For a bilayer we obtain valence bond 
phases with interlayer singlets, with alternating planar singlets, and two entangled spin-orbital 
(ESO) phases, in addition to the antiferromagnetic and ferromagnetic order. Possibility of such 
entangled phases in a monolayer is under investigation at present. 
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It has been shown that quantum fluctuations are en- 
hanced near the orbital degeneracy and could suppress 
long-range order in the Kugel-Khomskii (KK) model 
called below the d 9 model. This model was introduced 
long ago for a perovskite KCuF 3 , a strongly correlated 
system with a single hole within degenerate e g orbitals at 
each Cu 2+ ion. Kugcl and Khomskii showed that orbital 
order can be stabilized by a purely electronic supcrex- 
changc mechanism. This happens for strongly frustrated 
orbital superexchange Q , and columnar Ising-typc of or- 
der is obtained Q in the two-dimensional quantum com- 
pass model. This model exhibits nontrivial symmetry 
properties which may be employed to perform efficient 
calculations for square compass clusters 

Orbital order occurs in a number of compounds with 
active orbital degrees of freedom, where strong Coulomb 
interaction localizes electrons (or holes) and gives rise to 
spin-orbital superexchange When spin and orbital 
pscudospins couple to each other, their order is usually 
complementary — alternating orbital (AO) order accom- 
panies ferromagnetic (FM) spin order, and ferro-orbital 
(FO) order coexists with antifcrromagctic (AF) spin or- 
der. However, the above Goodenough-Kanamori rules, 
see also are not satisfied in cases when spin-orbital 
entanglement (SOE) dominates Q, as for instance in the 
spin-orbital d 1 model on a triangular lattice Q. 

The spin-orbital superexchange KK model for Cu 2+ 
(d 9 ) ions in KCuF 3 with S — 1/2 spins and e g orbitals 
descibed by r = 1/2 pscudospin was derived from the de- 
generate Hubbard Hamiltonian with hopping t, intraor- 
bital Coulomb interaction U and Hund's exchange Jh li- 
lt describes the Heisenberg SU(2) spin interactions cou- 
pled to the orbital problem by superexchange J = At 2 /U, 



n 



/.•-■Ml™ 



7 7 

T T- 

'» 3 



4 



in 



where {r%, ra, ^3} depend on r\ = Jh/U [9J, and 7 = a, b, c 
is the bond direction. In a bilayer two ab planes are 



connected by interlayer bonds along the c axis EH (a 
monolayer has only bonds within a single ab plane). Here 
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are projection operators on a triplet (singlet) configura- 
tion on a bond (ij), and t// are the orbital operators for 
bond direction 7 = a, 6, c. They are defined in terms of 
Pauli matrices {af , erf} as follows: 
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Finally, E z is the crystal-field splitting which favors ei- 



ther 



y 2 (if E x >Q)otz = 3z 2 - r 2 (if E z < 0) 



orbitals occupied by holes. Thus the model Eq. (1) de- 
pends on two parameters: E z /J and r\. 

The spin-orbital model Eq. (JTJ describes also Cu02 
planes in La2Cu04, where indeed U ^> t and large 
E z /Jh — 0.27 favors holes within x orbitals 0|. The su- 
perexchange between Cu 2+ ions ~ 0.127 cV reproduces 
there the experimental value. In this paper we consider 
the model Eq. (JTJ for K3CU2F7 bilayer compound where 
nearly degenerate e g orbitals are expected. It has been 
shown that the magnetic state of K3CU2F7 is described 
by interlayer valence bond (VB) phase stabilized by FOz 
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order with z orbitals occupied by holes 

We show below that the bilayer spin-orbital d 9 model 
Eq. {TJ describes a competition between different types 
of spin-orbital order. Consider first \E Z \ — > 00, where 
depending on the sign of the crystal field E z we get ei- 
ther FOz or FOa; configuration with (rf) = ±1/2 and 

(r° ) = Tl/4. After inserting these values into Eq. 
((T|) one finds the Heisenberg model describing cither an 
AF bilayer (E z — > —00) or two independent AF planes 
(E z — > 00) as in La2Cu04. In the limit of rj — > (l/3)~, 
the coefficient n = 1/(1 — 377) diverges and at large 
77 > 0.26 one finds fully FM configuration with AO order. 

The simplest approach is a single-site mean field (MF) 
approximation applied to the model Eq. ([T]). It excludes 
any spin fluctuations so the spin projectors n^*| (lit-,) 
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can be replaced by their mean values, where the depen- 
dence on the bond (ij) reduces to direction 7 in phases 
with translationally invariant magnetic order listed in Ta- 
ble |U the G-AF phase, the G-AF phase with AF planes 
and FM interplane bonds, the A-AF phase with FM 
planes and AF interplane bonds and the FM phase. 

In the orbital sector we apply then the MF decoupling 
for the products {r^r^ } along the axis 7: 



T 7 7 



+ t7(^ 7 )-(t7)(t7 
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As order parameters we take t a = (rf) and t c = (rf ) for 
a chosen site i = 1 (which is sufficient in orbital sector 
as t b = —t a — t c ) and we assume two orbital sublattices: 
each neighbor of the site i is rotated by ir/2 in the ab 
plane meaning that (r^j) = t b ^ a \ The self-consistency 
equations can be solved analytically (see Ref. [l(J) and 
the phase diagram of Fig. [lja) is obtained by compar- 
ing the ground state energies for different points in the 
(E z /J,rj) plane. One finds two classes of solutions: (i) 
uniform orbital configurations (t c = ±1/2, t a ^ = Tl/4) 
for global FO order, and (ii) nontrivial AO order with 
orbitals staggering from site to site in ab planes. 

For rj = we have only two AF phases, see Fig. QJa): 
G-AFz for E z < -J/4 and G-AFx for E z > - J/4, with 
different FO orders involving z or x orbitals, respectively. 
Because of the planar orbital configuration in the lat- 
ter G-AF phase one finds no interplane spin coupling 
and thus this phase is degenerate with the G-AF one. 
For higher rj the number of phases increases abruptly by 
three phases, all with AO configurations: the A-AF, G- 
AF/AO and G-AF/ AO phase. Surprisingly, the AO ver- 
sion of the G-AF phase is connected neither to FOz nor 
to FOx order in an antiferromagnet, excluding the mul- 
ticritical point at (E z /J,rj) — (—0.25,0), and disappears 
completely for 77 w 0.118. The G-AF/AO phase stays on 
top of uniform G(G)-AF phase, lifting their degeneracy 
at relatively large 77 and then gets replaced by the FM 
phase which always coexists with AO order, so one can 
conclude that the G-AF/G-AF degeneracy is most easily 
lifted by turning on the orbital alternation. On the op- 
posite side (for E z < 0), the G-AFz phase is completely 
surrounded by A-AF phase with AO order. In the A- 
AF phase the AF correlations in the c direction survive 
despite the overall FM tendency when 77 grows. This fol- 
lows from the orbitals' elongation in the c direction which 



TABLE I. Mean values of triplet and singlet spin projection 
operators for a bond (ij) in the ab plane and along the 
axis c in magnetic phases with long range order, see Fig. Q] 
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FIG. 1. Phase diagrams of the d 9 bilayer model (1) obtained 
in: (a) a single-site MF, and (b) a cluster MF. In (a) shaded 
(green) area indicates phases with nontrivial AO order while 
the other ones have FO order. In (b) light shaded (yellow) 
area marks singlet phases with spin disorder and dark (or- 
ange) shading indicates phases with SO entanglement. 



stabilizes interplane singlets in a better cluster MF ap- 
proach, sec below. Finally, the FM phase is favored for 
any E z if only r] is sufficiently close to 1/3, as expected. 

In a better cluster MF (or Bethe-Peierls- Weiss) ap- 
proach, introduced to capture the effects of quantum fluc- 
tuations, one divides the bilayer square lattice into sep- 
arate cubes containing 8 sites each and treats the bonds 
inside a cube exactly, and the bonds connecting different 
cubes in MF. This approach has at least three advan- 
tages over the single-site MF: {%) spins can fluctuate, (ii) 
elementary cell can double, and (Hi) we can have inde- 
pendent spin-orbital order parameter. The MF leads in 
a cluster to three order parameters: magnetic (s) = Sf, 
orbital t a( - b \ and on-site SOE r a ^ = (Sfrf 1 ^) - st a ^ . 

The self-consistency equations take rather complicated 
form (for details see Ref. [l(|) and can be solved only 
numerically by time-consuming iterative Lanczos diago- 
nalization of a cluster combined with updating the MFs. 
In orbital sector apart from the AO order described ear- 
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FIG. 2. Order parameters spin s, orbital ' and on-site 
SOE r a ' 6 ' for r\ — 0.15 and for increasing: (a) _B 2 in a bilayer 
in the ESO, EPVB and PVB phase, for T = 0; (b) T in 
a monolayer in the ESO, AO-PM and PM phase, for E z = 
-0.6J. 



lier, we consider configurations where orbitals within a 
cluster break the symmetry between a and b directions 
but the neighboring clusters are rotated by tt/2 in the 
ab plane, so globally the symmetry is preserved and the 
elementary cell is doubled. In the spin sector we consider 
the same configurations as in a single-site approach. 

Including spin fluctuations in the cluster MF approach 
stabilizes the G-AF phase with x orbitals over the FM 
one (E z > 0) but suppresses it when z orbitals are filled 
by holes (E z < —J/4), and gives instead three singlet 
VB phases called: PVB (plaquette VB), VBz and VBm, 
see Fig. H{b). VBz phase replaces G-AFz phase shown in 
Fig. QJa) and involves interplane singlets accompanied by 
FOz configuration. This phase was observed in K3CU2F7 
by Manaka et al. [ll| — here we explain it for realistic 
r] ~ 0.14. The VBm phase is very similar to VBz but 
with slightly modified FO order by an AO component 
increasing toward the A-AF phase. Transition from VBm 
to VBz is of the second order. In the PVB phase spin 
singlets are pointing uniformly in a or b direction within 
the cluster and the elementary cell is doubled. 

A different class of phases involves SOE — these are 
the ESO, EPVB and PVB-AF phase. All of them exhibit 
SOE but only the ESO and EPVB ones lie in the highly 
frustrated part of the phase diagram and have large on- 
site entanglement r a ^ b \ as shown in Fig. [Ha). The 
PVB-AF phase connects PVB and G-AF phases by sec- 
ond order phase transitions and is characterized by fast 
changes in orbital order and appearance of global mag- 
netization. The ESO phase has no magnetization and 
FO order is here much weaker than in the VBz phase. 
When E z grows, the ESO phase does change continu- 
ously into the EPVB configuration, being an entangled 
precursor of the PVB phase, with doubling of the unit cell 
and finite AF order which vanishes smoothly approach- 



ing the PVB phase. Additional calculations described in 
[lOj show that these entangled phases are absent if one 
assumes that (SiT^^) factorizes, i.e., r a< ^ = 0. 

Using the same cluster MF approach as above one can 
easily study the phase diagram of the KK model for a sin- 
gle layer at finite temperature T. At T = one finds the 
AF, FM, and PVB phases together with an ESO phase 
between the AF and FM phases. Turning on the ther- 
mal fluctuations we have found that typically the orbital 
order is much more robust than the magnetic one and 
the orbital configuration compatible with lattice geom- 
etry can greatly stabilize spin order. In Fig. [2jb) we 
present the thermal evolution of the order parameters 
{s,t a{ - b) } and on-site SOE parameter r a ^ in the ESO 
phase which melts and ends up as an ordinary paramag- 
netic (PM) phase. More details and the phase diagrams 
will be reported elsewhere. 

Summarizing, we have shown that spin-orbital entan- 
glement leads to exotic types of order which are stabi- 
lized by quantum fluctuations both in bilayer and mono- 
layer systems.. They emerge from highly frustrated 
spin-orbital supcrexchange and could be discovered only 
within a cluster mean field approach. 
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